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Using N. Euler's theorem on the integrabihty of the general anharmonic oscillator equation [i|], 

Z^ ■ we present three distinct classes of general solutions of the highly nonlinear second order ordinary 

, I differential equation ^^ + /i (t) ^ + /2 (t) x + fs (f ) x" = 0. The first exact solution is obtained from 

O ,■ a particular solution of the point transformed equation d^X/dT^ + X" (T) = 0, n ^ {~3, —1, 0, 1}, 

.^^ ' which is equivalent to the anharmonic oscillator equation if the coefficients fi{t), i = f,2,3 satisfy 

, an integrability condition. The integrability condition can be formulated as a Riccati equation for 

CO ■ /i (i) and 7-775- ^ respectively. By reducing the integrability condition to a Bernoulli type equation, 

two exact classes of solutions of the anharmonic oscillator equation are obtained. 

-(— > 



INTRODUCTION 



The anharmonic oscillator is a physical system generalizing the simple linear harmonic oscillator ^^ +ujlx{t) = 0, 
where x{t) is the position coordinate, t is the time, and ujq is the oscillation frequency. In general, the time evolution 
of the space variable x of the anharmonic oscillator is governed by the following nonlinear second order differential 
^ ■ equation [H, ^ 

^*^ I H T cIt 

^. ^+/l(t)^+/2Wx + /3(i)x"=/4(t), (1) 

'^nT . where /j (i), i = 1, 2, 3, 4, and x are arbitrary real functions of t defined on a real interval / C JR, with fi it) and x (t) 
^J ■ g C^^). The factors fi (t) are physically interpreted as follows: /i (i) is a damping factor; /2 it) is a time dependent 

^ , oscillation frequency coefficient; /a (t) is the simplest possible anharmonic term; f^ (t) is a forcing term, and n is a real 

constant [3|- The equation of motion of the anharmonic oscillator is strongly nonlinear, and when the anharmonicity 
term /a [t) x" is small, its solutions can be obtained by using perturbation theory. If the anharmonicity is large, then 
other numerical techniques need to be implemented. 

The anharmonic oscillator equation Eq. ([1} with specific values of the exponent n can be used to model many 
Cd ■ different physical systems. For n = 3/2, one obtains the Thomas and Fermi atomic model [3, Hli while the case 
n ~ —3 corresponds to the Ermakov Q, or Pinney [a], equation. For n = —1 one obtains the Brillouin electron beam 
focusing system equation [3, Q , and n = 3 gives the case of the Duffing oscillator [0] . 

An interesting particular case of the general anharmonic oscillator equation Eq. ([1]) is the Ermakov-Pinncy equation 
(EPE), which is a well-known exainple of a nonlinear second order differential equation with important physical 
applications (we refer the reader to [l3, \lM foi' ^^ historical development and an excellent review of the properties of 
the EPE equation). The EPE is endowed with a wide range of physic al ap plications, including quantum cosmology 
[13 | , dynamics of scalar field cosmologies and the braneworld scenario [13| , quantum field theory [IJ, [l^ , nonlinear 
elasticity [16[, nonlinear optics [Ojllll: description of the wavefunction of Bose-Einstein condensates (BEC) at the 
mean-field level [1^, the envelope of the electric field in nonlinear optics [20|, amongst others. In this context, the 
EPE provides an effective description for the time-dependence of the relevant spatially dependent field, typically 
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being associated with its width both in the EEC [2l|, [22] and in optical settings [23| . The mathematical analysis and 
structure of the EPE have been extensively discussed in |24| - [27| . 

Note that for generic values of the coefficients, Eq. ([1]) is equivalent to a third order autonomous dynamical 

system, which generically admits no closed form general solution [4I. The mathematical properties and applications of 

articular forms of Eq. ([T]) have been widely investigated, such as, the partial integrability of the anharmonic oscillator 

, the time-dependent driven anharmonic oscillator and its adiabaticity properties |28| , toroidal p-branes, anharmonic 
oscillators and (hyper) elliptic solutions [29[, conformal mappings and other power series methods for solving ordinary 
differential equations [30| , and the anharmonic oscillator in the context of the optimized basis expansion [3l| . The 
Painleve analysis of Eq. ([1]) was performed in |32| . Specific transformation properties of the anharmonic oscillator 
were considered in [l| , where an excellent review of the Lie symmetries approach to Eq. ([3]) can also be found. 

The most general conditions on the functions /i, /2 and /s J[or which Eq. ([3]) may be integrable, as well as conditions 
for the existence of Lie point symmetries, were obtained in [l|. Time-dependent first integrals were also constructed. 
The main results of [l| are that if n ^ {~3, —1, 0, 1}, then Eq. ([3]) can be point transformed to an equation of the 
form cPX/dT'^ + X" (T) = 0, can be linearized as tPX/dT^ -|- fcj = 0, fc2 € ^\0, and it admits a two-dimensional Lie 
point symmetry algebra. 

It is the purpose of the present paper to obtain, by using the results of [l| , some classes of exact solutions of the 
anharmonic oscillator Eq. ([T]) without the forcing term. The first solution is obtained by considering a particular 
solution of the point transformed equation (PX/dT^ + X" (T) = 0, equivalent to the initial anharmonic oscillator 
equation. The integrability condition obtained in [l| for the anharmonic oscillator can be formulated in terms of a 
Riccati equation for the fi{t) and -rj7)-if terms, respectively. By imposing some specific constraints on the coefficients 
of the Riccati equation, namely, by requiring that the Riccati equation can be reduced to a Bernoulli equation, two 
distinct classes of exact solutions of the anharmonic oscillator equation with zero forcing term are obtained. In the 
analysis outlined below, we shall use the generalized Sundman transformations X (T) = F (t, x) and dT = G (t, x) dt 
[33l435| . The latter have been widely applied in the literature |33l . |34{ . namely, in the study of the mathematical 
properties of the second order differential equations, and the third order differential equation 



d T d X dx 

^ + /ii (0 ^ + h2 (i) -^ + ^'3{t)x + hi (t) = 0. 



(2) 



The present paper is organized as follows. Three distinct classes of general solutions of Eq. ([T]) without the forcing 
term, which explicitly depict the time evolution of the anharmonic oscillator, are presented in Section |lll We discuss 
and conclude our results in Section IlIII 

II. EXACT INTEGRABILITY CASES FOR THE ANHARMONIC OSCILLATOR 

In the present Section, by starting from the integrability condition of the anharmonic oscillator equation obtained 
in [l|, we obtain three cases of exact integrability of the anharmonic oscillator without forcing. 

A. The integrability condition for the anharmonic oscillator 

In the following we assume that the forcing term f4{t) vanishes in Eq. Q. Hence the latter takes the form 

^ + fiit)^ + f2it)x + hit)x-^0. (3) 

An integrability condition of Eq. ^ can be formulated as the following: 

Theorem [l|. If and only if n ^ {—3, —1, 0, 1}, and the coefficients of Eq. ([3]) satisfy the differential condition 



/2W 
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1 dVs 



n-t-4 



1 dfs 



h{t) dt 



n~ 1 

(n + 3)' 



1 dh 

hit) dt 



flit) 



+ 3 hit) dt^ („ + 3) 
with the help of the pair of transformations 

X{T) - Cxit)ff^'it)e^^^'f 
T{x,t) - C^ f /3^(0e(^)^'^i('^)'^'^de, 



■^^-^^/?W. (4) 



n + 3 dt (n + 3)^ 



(5) 
(6) 



where C is an arbitrary constant. Eq. ([3]) can be point transformed into the second order differential equation for 

X[T), 



<PX 



The general solution of Eq. ([7]) is given by 

T = To + e 



X" (T) == 0. 



dX 



2 (Co ^ 



.+1 ) 



n^-1, 



(7) 



(8) 



where Tg and Co are arbitrary constants of integration. For convenience, we have denoted T = T±, Cq ~ Co± 
To = To± and e = ±. 

By substituting the integrability condition given by Eq. (JU into Eq. ^, we obtain the following integrable differ- 
ential equation 



1 1 d^fs 



71 + 4 



d X . . dx 



1 dh 
/3(i) dt 



n — 1 



1 rf/3 



(n + ar L/3(i) rfij 



/l(0 



— — ^^^^^/f(<) \x + h{t)x-=0,n^ {-3,-1,0,1}- (9) 



+ 3 rff (n + 3)' 



B. A particular exact solution for the anharmonic oscillator equation 

The general solution of Eq. ([7]) can be given as 

' n + 3 n + 2 X"+i 



™^-l/^^^^^T^i.F- 



1,- 



2(n+l)'n+l' Co(n + l) 
where 2-^1(0, b; c; d) is the hypergeometric function. A particular solution of Eq. ([7]) is given by 



n^-1. 



X{T) = [e{T-To)]-- 



2(71+1) 



(10) 



(11) 



where we have defined X (T) = X± (T), and we have taken the arbitrary integration constant as zero, Cq = 0. In 
order to have a real value of the displacement x{t) one must impose the condition 7i < — 1 on the anharmonicity 
exponent n. From Eqs. (O and ([TT|) . we obtain the result 



x{t) = -[e{T-To)]- 



2(71+1) 



f3"^''{t)e- — 



(12) 



where we have denoted x (t) = x± (t), for simplicity. 

By inserting Eq. (|6]) into Eq. (|T2|) yields the general solution of Eq. (j9]) describing the time evolution of anharmonic 
oscillator. Therefore we have obtained the following: 

Corollary 1. The anharmonic oscillator equation Eq. © has the particular solution 



X (t) = xo 



C"^ / /3^(C)e(^)^'^^(^)''*de-To 



/a -^' {t)e-^^'f'^^')''^,n<-l. 



(13) 



where we have defined xq = C ^ 



2(n+l) 



C. Second integrability case for the anharmonic oscillator equation 

Now, by rearranging the terms of Eq. (|4]) yields the foUowing Riccati equation for /i (i) given by 

^=ait)+bit)h{t)+cit)f!{t), 



where the coefficients a{t), b{t) and c{t) are defined as 

3 + n_., 1 d^fs 



a{t) 
bit) 
c{t) 



2 
1-n 


2(3 + n) 

f + 71 
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4 + n 
2/3(0 dt^ ' 2 (3 + n) 
1 #3' 



1 rf/3 



1 2 



h it) dt 



h it) dt 



3 + n 



We consider now that the coefficient a it) of Eq. (|T4)) vanishes, so that Eq. (jTS)) can be written as 



/2(i) 



1 1 d\h 



4 + n 



3 + n/3(f) dt2 (3 + „)2 



1 d/3 



/3 (t) dt 



Hence the Riccati Eq. p4|) becomes a Bernoulh type equation, with the general solution given by 

hit) 



f!^ it) 



<-i + 3+;r J h 



(^)d</. 



where Ci is an arbitrary constant of integration. 

By substituting Eqs. (fT8|) and (|f9p into Eq. ([3|), the latter yields the following differential equation 



di2" 



/J^^(i) 



c^i+^r/3'^(^)' 
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4 + n 



3 + n/3(i) dt2 (3 + nf 



1 d/3 



h it) dt 



(14) 

(15) 

(16) 
(17) 

(18) 



(19) 



X + /3 it) x" = 0. (20) 



by 



Therefore we have obtained the following: 

Corollary 2. The general solution of Eq. ([20)1 . describing the time evolution of the anharmonic oscillator, is given 



xit) = .To 



c— 



fz^'iOe 



' {* 



'^i + ^S*S: 



2(3 + »i) 



{^)dii 



d^-To 



x/3"^Mt)e 



l + n ,-,(, ,2(3+71 



^{-3,-f,0,f}. 



(21) 



D. Third integrability case for the anharmonic oscillator equation 

Now, by introducing a new function u it) defined as 



it) 



1 dh 



or, equivalently, 



/3 it) dt ' 

/3(t) = /o3eJ■'"(^)'^^ 



(22) 



(23) 



where /03 is an arbitrary constant of integration, after substituting u {t) into Eq. ([4]) yields the following Riccati 
equation for u (t), given by 



^ = ai (t) + h {t) u (t) + ci (t) u^ {t) 



where the coefficients are defined as 






6 + n at 



ci(t) 



3 + n 
1 

3 + n' 



We consider that the coefficient ai (t) of Eq. ([M)) vanishes, as before, so that Eq. (^5]) can be written as 

2(l + n) , 2 d/i 

•'2 (,'j - 77 ^/l [t + — — --77. 

(3 + n) (3 + ?i) di 

Then the Riccati Eq. (|24p becomes a Bernoulli type equation, with the general solution given by 



uit) = 



e3 






C-Jhf^^'''' 



where C2 is an arbitrary constant of integration. Thus, the coefficient /a (t) of Eq. ^ is readily given by 



r 






dV 



/3 (i) = /ose 
By substituting Eqs. I\28\i and (pO)) into Eq. ([3]), the latter yields the following differential equation 
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by 



Therefore we have obtained the following: 

Corollary 3. The general solution of Eq. (|3T|) . describing the time evolution of an anharmonic oscillator, is given 



X {t) = Xo/o, '+" 



^03 
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\/^ 



■ S^' fl{4-)d4> 
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_|qr^/P/l(*)<i*^ 



-+2/l(»A) 
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(32) 



III. CONCLUSIONS 



In the limit of a small function X (T), and by assuming that the constant n is large, n — > +00, in view of Eq. ([7]) 
we obtain a linear relation between X (T) and T (t), given by 



X{T)^e^2CoiT^To). 



(33) 



With the help of Ec^s. (O and ([6]), the approximate solution of Eq. (j3]), describing the time evolution of anharmonic 
oscillator is given by 



'{t) 



eV2C^ 
C 



c^ I /3"+-MOe^^)^'-^'*^'''^dC-ro 



/3-'-^(t)e-- 



■ f fl{<t>}d<P 



(34) 



With this approximate solution, once the functions /i (t) and /s (t) arc given, one can study the time evolution of the 
anharmonic oscillator for small X (T) and for a very large anharmonicity exponent n. 

in the present paper, by extending the results of [i[, where the first integral of Eq. ([3]) was obtained, we have 
obtained three classes of exact general solutions of Eq. ^, by explicitly showing that the Theorem obtained in [l| 
is very useful for obtaining the explicit general solutions of the anharmonic oscillator type second order differential 
equations. 

In order to have real solutions, the general solutions Eqs. ([T5|. (PT|) and (|5^ of the second order differential Eqs. 
(O; (HOI) and ([?!]) . respectively, must obey the condition n < — 1, thus leading to an anharmonic term of the form 
/3(t)/a;", n > 0. Such a term may be singular at x = 0. Note that in [301, the author has used the degree theory to 
remove a technical assumption in the non-resonance results in [37[ and to obtain a complete set of the non-resonance 
conditions for differential equations with repulsive singularities. In doing so, a nice relation between the Hill's equation 
^jf + T] {t) X (t) = and the EPE was established. This relation itself is useful in studying the stability of periodic 
solutions of Lagrangian systems of degree of freedom of 3/2. 

It is well-known that the second order ordinary differential equations with anharmonic term of the form f^{t)/x"' 
entail many problems in the applied sciences. Some examples are the Brillouin focusing system, and the motion of an 
atom near a charged wire. The Brillouin focusing system can be described by the second order differential equation 

(fix P 

— +a{l + cost)x{t)^——-, (35) 

where a and /3 are positive constants. In the context of electronics, this differential equation governs the motion of a 
magnetically focused axially symmetric electron beam under the influence of a Brillouin flow, as shown in [7| . From 
the mathematical point of view, this differential equation is a singular perturbation of a Mathieu equation 

—^ + (a-2qcos2t)x(t) = 0, (36) 

at'' 

where a and q are arbitrary constants. Existence and uniqueness of elliptic periodic solutions of the Brillouin electron 
beam focusing system has been discussed in [8|. Hence, the results obtained in the present paper could open the 
possibility of obtaining some exact solutions of non- linear differential equations of scientific or technological interest. 
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